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Abstract 

We have discussed dynamical behavior of strange quark matter components, in particular the 
effects of density dependent quark mass on the equation of state of strange quark matter. Dy- 
namical masses of quarks have been computed within Nambu-Jona-Lasinio (NJL) model, then we 
have done the strange quark matter calculations employing the MIT bag model with these dy- 
namical masses. For the sake of comparing dynamical mass interaction with QCD quark-quark 
interaction, we have considered the one-gluon-exchange term as the effective interaction between 
quarks for MIT bag model. Our dynamical approach illustrates an improvement for the obtained 
values of equation of state. We have also investigated the structure of strange quark star using 
Tolman-Oppenheimer-Volkoff (TOV) equations for all applied models. Our results show that the 
dynamical mass interaction leads to lower values for the gravitational mass. 
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I. INTRODUCTION 



Strange quark stars (SQS) are the most compact objects with a surface density p ~ 
10 15 ^_-, -^jch i s about fourteen orders of magnitude greater than the surface density of 



neutr on sta rs, while their c entra 
et al. 



density c ould be up to five times higher than that (Haensel 
20QZ1; Glendenning |2000fl ; Weber [19991 ] V It was first Itoh (1970) that, even before 
QCD full development, proposed SQSs which is made of strange quark matter (SQM). 
Later, Bodmer (1971) discussed the fate of an astronomical object collapsing to such a state 
of matter. 

The quark deconfinement hypothesis is one of the exciting steps in investigation for 
the building blocks of matter. Soon after predictions of quarks in theories and successful 
laboratory observations, many hadronic models were developed to describe the probable 
quark matter proposed at high energy regimes. In the 1970s, after formulation of QCD, 
perturbative calculations of the equations of state of SQM got form, but the ar ea of validity 
for these calculations was restricted to very high densities (Collins & Perry 19751 ] ) . The 
existence of SQSs was also discussed by Witten (1984), who conjectured that a first order 
QCD phase transition in the early universe could concentrate most of the quark excess in 
dense quark nuggets. Witten proposed that SQM composed of light quarks is more stable 
than nuclei, therefore SQM can be considered as the ground state of matter. 

An SQS would be the bulk SQM phase consisting of almost equal numbers of up, down 
and strange quarks, plus a small number of electrons to ensure charge neutrality. A typical 
electron fraction is less than 10 -3 and decreases fro m th e surface to t he ce nter of an SQS 



(Haensel et al. [20071 ]: Glendenning [20001 ]: Weber [19991 ]: Camenzind [20071]). SQM would 
have a lower charge-to-baryon ratio compared to t he nuc lear matter and can s how itself in 
the form of an SQS (Witten [198J]; Alcock et al. [1986) : Haensel et al. [1986| : Kettner et 



al. 



19951 ] ). The collapse of a massive star could lead to the formation of an SQS. An SQS 



may al so be formed from a neutron star and is denser than the neutron star (Bhattacharyya 
et al. [20061 ] ) . If sufficient additional matter is added to an SQS, it will collapse into a 
black hole. Neutron stars with masses of 1.5 — 1.8M with rapid spins are theoretically the 
best candidates for conversion to an SQS. An extrapolation based on this indicates that up 
to two quark-novae occur in the observable universe each day. In addition, recent Chandra 
observations indicate that objects RX J185635 — 3754 and 3C58 may contain SQSs (Prakash 
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et al. 



2003|). 0ther i nvesti gations also show that the object SWIFT J1749.4 - 2807 may 



be an SQS (Yu & Xu |2010|). 

The strange quark star, founded from quark matter theory consists of too many unsolved 
puzzles which are usually involved in physics of these relativistic objects. System complexity 
of these stars prohibit us from considering all physical and astrophysical properties simul- 
taneously, and it is possible that some parameters entering the equation of state do not 
represent specific physical properties. For example, in MIT bag model, one of the models 
used in this paper, when the researchers try to find an d fit the bag const ant ac cording to 
infor mation s achieved from big col liders (Jin Sz J enning 1997| ; Alford et al. 1998|; Blaschke 



et al. 



19991 ]; Burgio et al. [2002bj |; Begun et al. [20111 ]). we should keep this principle as a 



matter of fact that different parameters like temperature, electromagnetic intensity, density 
etc. are important enough on final interpretation for theoretical calculated bag constant. In 
this point of view even constant values of bag pressure no more can be considered purely as 
the energy density difference between the perturbative vacuum and the true vacuum. The 
role of bag constant for confining quark matter in comparison with gravity confinement for 
neutron matter may require more attention when we consider it for compact stars. There- 
fore it is better to consider the dynamical properties of the parameters for investigating of 
the properties of quarks. Many works have been done to adapt the theory of bag Model 
on physics of ultra-dense matter like using a density dependent bag constant (Burgio et 



al. 



2002a | ), u tilizing different values of coupling constants for one gluon exchange (Farhi 



& Jaffe [19841 ] : Berger & Jaffe [19871 ] ), o r inste ad consideri ng dyn amical mass as effective 
interaction between particles (Peng et al. jl999 |: Shao et al. 2011 1 ) . 

From perturbative QCD, we know that quarks at ultra high densities asymptotically 
interact. One way of considering the interaction is to assume that quarks exchange one gluon. 
Therefor we can add a term to equation of state that is characterized by a coupling constant. 
But constant values of this parameter will weaken the power of interaction in lower densities 
while in higher densities, it increases it. One method to solve the problem is to assume a 
density dependent quark ma ss as the effective interacti on. T his approach wa s inve st igated 
in references (Fowler e t al. 19811 ] ; Chakrabarty et al. 19891 ]; Chakrabarty 19911 ] . 19941 ] ; 
Benvenuto & Lugones [1995I ] ; Lugones & Benvenuto jl995| ). This has been done by adding 
a term to the rest mass which is characterized by a free parameter determined by stability 
conditions. They concluded that the density dependent mass is flavor independent and the 
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applied free parameter has the same meaning as the bag constant. Then by selecting one 
value of bag constant for all densit ies and flavors, they tried to obtain the equation of state 



of quark matter (Peng et al. 



19991 ]) . A better approach closer to current work, is t o find a 



solution for density dependent mass from Nambu-Jona-Lasinio (NJL) method (Carol 20091 ] ) . 
Carol calculated the equation of state and structure for hybrid stars within MIT bag model, 
while numerical values of density dependent mass entering in the energy equation had been 
obtained from dynamical calculations of mass in NJL model. These numerical values were 
entered directly in the pressure equation without considering density dependency. Quark 
masses and NJL constants were also approximate values. The bag constant in that work was 
density independent; therefore in addit ion to previous kno wn p roblems of const ant va lues 
for this parameter (Baldo et al. 120061 ]: Alford & Reddy |2003|; Alford et al. [2005| ). it 
misinterprets the meaning of the effective interaction in some densities. 

In Our previous w ork we considered a hot strange star just after the collapse of a super- 
nova (Bordbar et al. 2011al ]). at finite temperature with a density dependent bag constant. 
The calculations for the structure properties of the strange star at different temperature 
indicates that it's ma ximum mass decreases by increasing the temperature. In another work 
(Bordbar & Peivand 2011bj ). we concentrated on the calculation of a bulk of spin polar- 
ized SQM at zero temperature in the presence of a strong magnetic field. We computed 
structure properties of this system and found that the presence of a magnetic field leads to 
a more stable SQS when compared to the structure properties of an unpolarized SQS. In 
present paper, we investigate the quark matter equation of state and the strange quark star 



structure following Carroll (Carol 20091 ] ) . We base our calculations on MIT bag model, and 
after following NJL formalism we extrapolate a density dependent equation from numerical 
values of dynamical mass obtained using NJL method. In Sec. [Til the required equations 
for the MIT bag model are written, the same has been done for NJL model. In Sec. Ill CI 
we describe the formalism applied in this article, and after solving TOV equations in Sec. 
IHIl we calculates SQS structure for our method. 



II. CALCULATION OF EQUATION OF STATE FOR SQM 



In this section, we calculate the equation of state of strange quark matter (SQM) using 
MIT and NJL methods as well as MIT method with the dynamical mass. At first, we 
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introduce these three models in three separate sections, then we give our results for the 
energy and the equation of state of SQM in sec. Ill Dl 



The MIT Bag Model 



Total energy of a bulk of dec onfine d up (u), down (d) a nd stran ge (s) quarks withi n MIT 



bag model is as 
Berger & Jaffe 



ollows (Witten [198 



19871 ] : Glendenning 



; Farhi & Jaffe |l984j : Baym 



19851]; Baymet al. [1985j ; 



19901 ]: Maruyama et al. |2007j ): 

e = e u + e d + e s + B. 



(1) 



In Eq. (Q]), B is the bag constant, and 



£ f (Pf) 



3m f 4 
8vr 2 



ni 



f 



Xf (2x 2 f + lj ^\Jl + x 2 ^j — arcsinh Xf 

3 r / c o 

1 + x 2 



(2) 



xj — - \xj y^j i + XJ) — arcsinh Xf 

where / denotes the flavor of the relevant quark, a c is QCD coupling constant and the 
following term demonstrates the one-gluon-exchange interaction. In above equation, Xf is 
defined as follows, 

Xf = k F (fi/nrif, (3) 
where the Fermi momentum &p W is given by 

2\V3 



PfK 



(4) 



For the ba g constant (B ), we use a density dependent Gaussian parametrization (Burgio et 



al. 



2002a) : Baldo et al. |2006|): 



B (p) = ^00 + (B - B^) exp[-/3 (p/p ) 2 ] (5) 

with B^ = B (p = oo) = 8.99 MeV/fm 3 , B = B (p = 0) = 400 MeV/fm 3 and /3 = 0.17. In 
SQM, the beta-equilibrium and charge neutrality conditions lead to the following relation 
for the number density of quarks, 



P = Pu = Pd = Ps 



(6) 



From the total energy, we can obtain the equation of state of SQM using the following 
relation, 
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B. The Nambu-Jona-Lasinio Model 



Here we give a brief introduction regarding the calculations in the Nambu-Jona-Lasinio 
(NJL) method. For NJL model, we use a common three flavor lagrangian adopted from 
(Rehberg et al. 19961 ]) which preserves chiral symmetry of QCD, 



C = q(iYdf, -m )q + GJ2 [(^kqf + (qilshqY 



(8) 



fc=0 



K [det f (q (1 + 75) q) + det f (q (1 - 75) <?)] ■ 



In adopted lagrangian, q denotes quark field with three flavors u, d and s, and three colors. 
m = diag(rriQ, rag, tuq) is a 3 x 3 matrix in flavor space. And ( < k < 8 ) are the 
U(3) flavor matrices. We restrict ourselves to the isospinsymmetric case, 



have pick ed up the parameters from references (Kunihiro 19891 ] ; Ruivo et al. 1999| ; Buballa 



o- 



We 



& Oertel 



1999]) which are fitted to the pion mass, the pion decay constant, the kaon mass 



and the quark condensates. 

NJL model is an unrenormalizable method with divergent integrations. To prevent the 
divergence, we need to introduce some breaking points for upper limit of integrals which 
satisfy the physics ranges of our problem. It is usually done by choosing a proper cut-off. 
In present paper, the adopted cut-off is named Ultra-violet cut-off that indicates restoring 
of chiral symmetry breaking, A = 602.3 MeV. G and K are coupling strengths that read, 
GA 2 = 1.835, KA 5 = 12.36. The rest mass of s quark is m s Q = 140.7 MeV, and there is 



m 



m 



5.5 MeV for u and d quarks. The baryon number density is given by 



1 



Pb 



1 



-n B 



n u + n d + n s ) 



(9) 



3 3 

where rij = (q\qi)- Within mean field approximation, the dynamical mass is calculated by 
the following gap equation, 



mi = m\- AG (^^) + 2K (q^qj) (q k q k ) 



(10) 



In the above equation, we need to calculate permutation of all quark flavors. The quark 
condensate in Eq. (fT0l) reads 

3 r A 



P Z d P 



mi 



mf + p 2 



(11) 
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and Pfi, Fermi momentum of quark i, is obtained from the following relation, 

P Fi =U 2 ni)K (12) 



Equations ( TTOj) and (iTTj) have self consistent solutions. It means that for a given number 
density, n.j, we should calculate quark condensate and substituting the corresponding value in 
Eq. ( fTOl) to reach a consistent result of the dynamical mass after doing the iteration process. 
In Fig. CE} we have plotted the results of density dependent mass for u, d and s quarks as 
a function of density. As it is clear from Fig. [H quark masses vary from current masses 
(5.5MeV for u and d quarks, and 140.7MeV for s quark) at high densities to constituent 
mass at near zero densities (368.7 MeV for u and d quarks, and 550MeV for s quark). 

The solution via mean field approximation forces us to stabilize equations by diminishing 
energy density and pressure in vacuum. This is satisfied by entering a pa ramet er which has 



the same meaning of bag constant in MIT bag model (Buballa & Oertel [1999]): 

'3 rA 



B = E (-2/ p 2(1 p Up 2 + m "- "Jp 2 + mi l) - 2G fe) 2 ) 

i=u,d,s ® ) 

+ AK (uu) (dd) (ss) (13) 
Now we can calculate the equation of state of SQM in NJL model, 



p=-e+ E n^Ppj + ml (14) 

i=u,d,s 

where 

e= E 4 [ PFi p 2 dp^p 2 + mi -(B- B ) . (15) 



i=u,d,s 



Parameter B is the bag pressure, which is explained by Buballa (2005), and is a dynamical 
consequence of the mean field solution, not a parameter inserted by hand, as was done in MIT 
bag model. It is shown in Fig. [H matter in NJL method acquires dynamical mass in nonzero 
baryon densities, but in MIT bag model, the given mass remains constant for all densities. 
Consequently, this will lead to dissimilar chiral symmetry behavior as density changes. In 
NJL model, since quarks acquire dynamical mass, the chiral symmetry spontaneously breaks 
in lower densities, while in MIT bag model, it will happen physically when quarks change 
their directions by hitting the bag (what is not considered theoretically in ordinary MIT 
bag model). The bag constant versus density is presented in Fig. [2] for our used models. 
It is apparent from Fig. |2] that chiral symmetry in our calculations is fully restored in the 



densities greater than p ~ 2.5 fm~ 3 . It is also important to mention that vacuum in MIT 
bag model is totally free of particles (flow of particle's wave function is restricted by the 
confinement), while in NJL model no confinement is produced. In other word, the vacuum 
in NJL model is made of paired quasi-quarks that lower the energy density of particles in 
comparison to MIT bag model. From the above di scussio ns, it seems reasonable to add an 



effective bag constant to energy equation (Buballa [2005| ). 



Bo B | nii=nd=rajj= o, 

B eff = B-B . (16) 

From Fig. [2J it seems that the effective bag constant diminishes at zero density. Then the 
correct interpretation for the effective bag constant is the energy per volume needed to fully 
break quark-antiquark pairs in order to completely restore chiral symmetry at ultra high 
densities. Even the maximum value of dynamical NJL bag constant is smaller than that 
of MIT's one, because it reduc es the energy per particle due to quark-antiquark pairing at 



lower densities (Buballa [2005 | ) . Fig. |2] shows that the decreasing rate of MIT bag constant 
is higher than that of NJL. This indicates that MIT bag model does gross approximation 
over physics of matter in middle and higher densities (p > 0.8 /m~ 3 ). Therefore, the density 
dependent bag constant should be corrected by another higher density sensitive parameter. 
This could not be achieved by a one gluon exchange term that considers the interaction with 
a constant strength in all energy regimes. Fig. |2]indicates that at the density p ~ 0.45 fm~ 3 , 
there is a cross point for the effective bag constant of NJL model and the bag constant of 
MIT model. As it is mentioned in above discussions, the bag pressure is the energy needed to 
confine particles where effective bag constant is energy needed to destabilize quark-antiquark 
pairs. Now, we can suggest that the hadron-quark phase transition can takes place at th e 



density p ~ 0.45 fm 3 . This is in good agreement with the results of others (Heinz 2001]; 



Heinz & Jacob 



2000^ . 



C. MIT bag model with dynamical mass 

In MIT bag model with dynamical mass, we consider the effect of dynamical behavior of 
the quark mass in calculating the equation of state of SQM within MIT bag model using 
NJL numerical mass results. In fact, we use the dynamical masses (Fig. [1]) for u, d and s 
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quarks in Eq. (|2j) instead of their fixed values. 



D. Our results for the energy and equation of state of SQM 

To distinguish numerous outcomes, we present the results of our calculations in three 
following models; 

• Model 1: MIT model by a density dependent bag constant and one gluon-exchange 
(a c = 0, 0.16, 0.5) as effective interaction. 

• Model 2: NJL model. 

• Model 3: MIT bag model by a density dependent bag constant, dynamical mass and 
one gluon-exchange (a c = 0, 0.16, 0.5) as effective interaction. 

Our results for the energy of SQM versus density calculated with above models have 
been plotted in Fig. O We see that for both MIT based calculations (models 1 and 3), 
at lower densities (p < 0.5 fm~ 3 ), the energy of SQM sud denly increases as the density 



decreases. This shows the concept of confinement (Buballa 2005]). For these two models, 
we also see that the energy of SQM gets to a minimum, then increases with a small rate. 
Fig. |3] shows that for model 1 and model 3, the energies of different coupling constants 
are nearly identical for densities p < 0.5/m~ 3 . However, they have a substantial difference 
as the density increases. We can see that at lower densities (p < 0.7 fm~ 3 ), the results 
of model 3 is considerably different from those of model 1. While this difference becomes 
small as density increases, specially for lower values of coupling constant, due to asymptotic 
freedoms which is the simple MIT bag model without interaction. From Fig. [31 it is seen 
that the energy of SQM in model 2 (NJL model) has finite values even at low densities 
showing no confinement. We also see that the energy of SQM from model 3 with smaller 
values of coupling constant is lower than that of model 2 for p > 0.7 fm~ 3 indicating a 
more stable state of quark matter at these densities. However, at very high densities, the 
difference between the results of these two models becomes negligible. 

In Fig. HI our results for the pressure of SQM have been plotted versus density. It can 
be found that for MIT bag model, the higher values of coupling constant leads to the stiffer 
equation of state for SQM. Fig. H] shows that by considering a dynamical mass for the quarks 
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(density dependent mass) in MIT model, we get the lower values for the pressure of SQM. 
For a c = 0.0, we see that the result of model 3 for the equation of state of SQM is nearly 
identical with that of model 1. It can be seen that for p > 0.6 fm~ 3 , our results for the 
pressure of SQM calculated by NJL model are nearly identical with those of model 3 and 
model 1 for a c = 0.0, while at lower densities, there is a considerable difference between 
them. 

In order to investigate the quark matter stability, the energy of SQM versus pressure has 
been plotted in Fig. [5j It is clearly seen that at zero pressure, the MIT bag model with 
a c = leads to the lowest value for the energy of SQM (950MeV fm~ 3 ) compared to other 
models. This value is comparab 



e with the result for the binding energy per particle of Fe 
(930MeV fm~ 3 ) (Witten 1984J ). This indicates that among different models used in this 



work, MIT model with a c = shows the most stable state of SQM. 



III. CALCULATION OF STRANGE QUARK STAR STRUCTURE 



The gravitational mass (M) and radius (R) of compact stars are of special interest in 
astrophysics. In this section, we calculate the structural properties of a strange quark 
star for our three models. Using the equation of state of strange quark matter for the 
models applied in this work, we can obtain M and R by numerically integrating the general 



relativistic equations of hydrostatic equilibrium, the 
equations, which are as follows (Shapiro & Teukolsky 



blma n-Qppenheimer-Volkoff (TOV) 



1983|) 



& m . 2 / % 
— — = Airr e ir) , 

dr 



(17) 



dp Gm{r)e{r) ( p(r) \ ( 4irr 3 p(r)\ f 2Gm{r) 



-i 



dr 



e ir <r 



m (r) c 2 



c 2 r 



where e (r) is the energy density, G is the gravitational constant, and 



m{r) = / Anr 2 e(r)dr 
J o 



(18) 



(19) 



has the interpretation of the mass inside radius r. By selecting a central energy density 
e c , under the boundary conditions P(0) = P c and m(0)=0, we integrate the TOV equation 



outwards to a radius r = i?, at which P vanishes. 
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In Fig. El we have presented our results for the gravitational mass of SQS versus the 
central energy density. Fig. [6] shows that at low energy densities, the gravitational mass 
increases rapidly by increasing the energy density, and it finally reaches to a limiting value 
(maximum gravitational mass) at higher energy densities. It is seen that the increasing rate 
of mass for Model 3 with higher values of coupling constant is substantially higher than 
those of other models. Table II III summarizes maximum gravitational masses of different 
applied models and the corresponding radii. As it seen from Table IHH we can conclude 
that using dynamical mass in energy equation and equation of state of SQM reduces the 
calculated maximum mass. This is in a good agreem ent with many observational data 



obtained from low mass compact stars (Zhang [20071 ] ) . It is interesting that in spite of 
considering dynamical mass as the effective interaction in MIT bag model (model 3 with 
a c — 0), we find the smaller SQS maximum mass in comparison to MIT bag model (model 
1) even without interaction (a c = 0). As it is obvious from Table ITTTt for models 1 and 3, the 
calculated maximum mass increases as strong coupling constant increases. This behavior 
demonstrates that ultra massive SQS with masses greater than M = 1.05M© are stars 
which are composed of highly interacting strange quark matter. We note that some studies 
indicate that there exist a big uncert ainty about mass and radius of ultra massive stars with 



M > 1.9M Q (Lattimer & Prakash [2010] ). These studies showed that the observed data 
of mass and radius for these stars, which commonly belong to X-ray stars, were wrongly 
calculated and the calculations were revised to the smaller values for mass and radius. The 
best example is PSR J 0751 + 1807 pulsar that initially was supposed to h ave a mass of 



M = 2.2 ± O.2M but recently revised to M = 1.26M (Lattimer & Prakash |2010| ). 

We have also plotted the gravitational mass of SQS versus radius for our three models 
in Fig. [7] It is seen that for all models, the mass increases by increasing the radius, but 
with different increasing rates for different models. Fig. [7] shows that for a given value of 
radius, the dynamical model (model 3) gives the smaller mass with respect to that of MIT 
bag model (model 1); however, for a c = 0, it is close to the result of NJL model (model 2). 
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TABLE I: Maximum gravitational mass (M max ) and the corresponding radius (R) for different 
applied models. 



M max (M ) R(km) 



Model 1; a c 


= 


1.43 


7.61 


Model 1; a c = 


0.16 


1.73 


8.17 


Model 1; a c = 


= 0.5 


2.6 


10.6 


Model 2 




0.98 


5.59 


Model 3; a c 


= 


1.05 


6.03 


Model 3; a c = 


0.16 


1.65 


6.98 


Model 3; a c = 


= 0.5 


2.3 


8.69 
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FIG. 1: Density dependent mass (m) versus density (p) obtained from dynamical NJL model. 
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FIG. 2: Bag constant as a function of density for NJL and MIT models. 
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FIG. 3: The energy per baryon versus density for models 1 and 2 (a), and model 3 (b). 
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FIG. 4: Pressure as a function of density for models 1 and 2 (a), and model 3 (b). 




FIG. 5: Energy per particle versus pressure for models 1 and 2 (a), and model 3 (b). 
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FIG. 6: Gravitational mass (M) in unit of solar mass (M sun ) versus central energy density (e c ) for 
models 1 and 2 (a), and model 3 (b). 




FIG. 7: Gravitational mass (M) in unit of solar mass (M sun ) versus radius (R) for models 1 and 
2 (a), and model 3 (b). 



18 



